Spin and Orbital Characters of Excitations in Iron Arsenide Superconductors 

Revealed by Simulated Fe L-Edge RIXS 
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We theoretically examine the orbital excitations coupled to the spin degree of freedom in the parent 
state of the iron-arsenide superconductor, based on the calculation in a five-band itinerant model. 
The calculated Fe Z/3-edge resonant inelastic x-ray scattering (RIXS) spectra disclose the presence of 
spin-flip excitations involving several specific orbitals. Magnon excitations predominantly composed 
of a single orbital component can be seen in experiments, although its spectral weight is smaller 
than spin-flipped interorbital high-energy excitations. The detailed polarization and momentum 
dependence is also discussed with predictions for the experiments. 



After intensive study of iron arsenides motivated by 
the discovery of high-temperature superconductivity, it 
has been recognized that both spin and orbital degrees 
of freedom are the key to understanding the physics of 
iron arsenides. 

Spin excitations in the parent and superconduct- 
ing compounds have been studied by inelastic neu- 
tron scattering experiments* 1 - In the parent compounds, 
anisotropic spin waves have been observed below the anti- 
ferromagnetic (AFM) transition temperature T^^-~— and 
the anisotropic spin excitations persist just above T^A^- 
These behaviors are consistent with theoretical calcula- 
tions within the random-phase approximation (RPA) for 
a five-orbital Hubbard models 

In the AFM ordered phase with orthorhombic struc- 
ture, an electronic anisotropy in the FeAs plane has been 
reported from scanning tunneling microscopy* optical 
conductivity^^ and angle-resolved photoemission spec- 
troscopy (ARPES) , 11 ' 12 Resistivity measurements have 
revealed an in-plane anisotropy not only in the or- 
thorhombic phase but also in the tetragonal phase above 
the structural transition in BaFe2As2 - 13 ' This suggests 
a possible presence of a nematic order forerunning the 
structure transition and the magnetic ordering. The 
ARPES measurements above Tn have shown the lift of 
de generacy between d zx and dy Z orbitals,"" which indi- 
cates a strong influence of the orbital degree of freedom 
on the nematic order and the electronic anisotropy below 
Tn»^ In fact, the in-plane anisotropy of optical conduc- 
tivity below Tn has been explained by taking into account 
orbital characters of interband excitationsi 17 i 18 

Since both spin and orbital characters contribute to 
the electronic states of iron arsenides, it is desired to 
detect both the excitations at the same time in the energy 
and momentum spaces. For this purpose, we propose 
resonant inelastic x-ray scattering (RIXS) tuned for the 
Fe L edge as a probe. The L-edge RIXS consists of two 
processes: the x-ray absorption process and the emission 
process by way of an intermediate state with core holes. 
The x-ray absorption is accompanied by creation of an Fe 
2p core hole and an electron in the 3d orbitals, where the 



2p core hole state has a hybridized spin with the orbital 
angular momentum due to the spin-orbit coupling. When 
the electron relaxed through the x-ray emission is not 
that excited from 2p, the resulting orbital occupation of 
3d electrons varies from the initial. 



Note that the spin of the relaxed electron may be either 
up or down since the core hole spin is hybridized. The 
final state may thus involves a spin flipj^r— This means 
that we can investigate spin excitations in addition to 
orbital ones through the L-edge RIXS. Recently, single 
magnon excitation has been observed in the Cu L-edge 
RIXS for cuprate superconductor La 2 - I Sr x Cu04^ 3 - 



In this study, we theoretically investigate magnon exci- 
tations and orbital excitations coupled to the spin degrees 
of freedom together with spectral features of Fe L 3 -edge 
RIXS in iron arsenides. Our calculations are performed 
for a five-band Hubbard model by using the RPA! and 
a fast-collision approximation! 20 ! 22 In the AFM phase, 
we find that the magnon excitations predominantly com- 
posed of single orbital component appear in L-edge RIXS 
spectra with a weak intensity as compared with or- 
bital excitations — the orbital excitations lie just above 
those in contrast with the case of cuprates, where d- 
d excitations are well separated from the magnon ex- 
citations^ The dominant orbital excitations above the 
magnon mode are found to be accompanied by the spin- 
flip process, producing composite excitations of the cou- 
pled orbital-spin degrees of freedom. The origin of the 
excitations is attributed to the spin-flip particle-hole ex- 
citation from occupied to unoccupied states. We also 
predict the polarization and momentum dependence of 
the Fe L 3 -edge RIXS prior to forthcoming experiments. 



We calculate the ground state of a mean-field five-band 
Hubbard model written with the ordering vector Q = 
(j^,0) [Nq = 2 for AFM, and 1 for paramagnetic (PM) 
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where and z/ represent an electronic orbital, e M is the 
on-site energy, f (A; fi, v) is the transition energy between 
Fe sites distanced by A, is the creation operator of 

the 3d electron with wave vector k and spin a (quantized 
along the x axis), U and J are the intraorbital Coulomb 
interaction and the Hund coupling, respectively, and the 
pair hopping is set equal to J. We use the same parame- 
ter set as in Ref. S [/ = 1.2 eV and J = 0.22 eV, which 
yields the magnetic moment 0.8/^b- The distribution of 
the local magnetization among different orbitals is as fol- 
lows: 3z 2 — r 2 , 0.11/xb; zx, 0.13/xb; yz, 0.19^b; x 2 — y 2 , 
0.07/iB; xy, 0.29/xb- The spin-density-wave order param- 
eter is defined as (ni Q ^ a ) = J2k( d k+iQ^a d ^va) (N 
is the number of k points in the first Brillouin zone of 
the PM system, and I ^ 0). 

In terms of the quasiparticles diagonalizing Hmf, 
we consider particle-hole pair excitations on the mean- 
field ground state \g) with energy E g and write the 
Hamiltonian as if p h for the particle-hole system. We 
define the bare susceptibility as Xq 1S2 » ^(q, q', = 

/i T 

4o+q+iQvai d ko+'Q/"il5>. where s i =t, I, +i - for 
the spin pair (oi,^) = (t,t), (1,1), (4-,t)» (t,l), 
respectively, and S2 is associated with (0-2,(72) as well. 
We here set 77 = 0.01 eV and q = q'. 

We calculate the dynamical susceptibilities within the 
RPA 



At 






.11 yif 



Xo 





x UyU 




X* 

(2) 

where the product of susceptibility and interac- 
tion between quasiparticles is taken as a matrix 
product represented in the orbital basis such as 
[XoVxLa = EpVA'r' Xo v x>Vt> v'Xu> x , and the nonzero 

elements of the interaction matrix are V~ 
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FIG. 1. (Color online) Spin excitation spectra of transverse 
mode for AFM case, (a) the total spectrum is plotted. Spe- 
cific components corresponding to xtZ is abstracted: (b) 

ji — xy and v — xy, (c) fi — zx and v = xy, and (d) 



fi — 3z 



and v — xy. 
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tions are taken with respect to the four orbital indices. 

In this study, we use E r as the unit of energy to eval- 
uate the excitation energy under the presence of renor- 
malization effects: The larger the effects, the smaller the 
value of E r , setting E r = 1.0 eV represents the case of 
no renormalization effects included. Comparing our data 
of spin excitation with inelastic neutron scattering spec- 
tra,^ we guess that the value of E r is around 0.5 eV. 

The calculated spectrum of the spin-transverse mode 
is plotted for the AFM case in Fig. [Ha): The spec- 
trum around (0,0) is shown here for RIXS study, al- 
though that around (it, 0) shows a strong intensity. In the 
low-energy region (< 0.15EV), there is a sharp magnon 
branch stemming from (0,0). Orbital analysis reveals 
that the magnon excitation is primarily associated with 
the icy-orbital component xty~*y as shown in Fig. QTb). 

xy xy 

This implies that the xy-orbital plays a significant role 
in the local magnetic moment, consistent with the cal- 
culated local magnetic moment, to which the xy orbital 
component shows the largest (about 40%) contribution 
while the others contribute less than 25% each. 

Above 0.15-E r , there are flat excitation branches 
around (0,0): the 0.2E r mode and 0.3-EV mode. The or- 
bital character of the 0.2E r mode is a type of zx — > xy, 
which means that xt^v is dominant [Fig. [IT c)]. On the 

zx zx 

other hand, that of the 0.3EV mode is 3z 2 — r 2 — » xy 
[Fig-Eld)]. 

The spin-longitudinal and charge modes are plotted 
in Fig. O Both show excitations around 0.3_E r . These 
0.3-E,. modes possess 3z 2 — r 2 —5- xy character as well as 
the spin-transverse mode. 

The calculated density-of-states shows a peak in the 
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(a) spin-longitudinal 



(b) charge 




FIG. 2. (Color online) Excitation spectra of spin-longitudinal 
(a) and charge modes (b). 




FIG. 3. (Color online) Situation setting. Here, x-ray incoming 
with momentum k by angle 6 is scattered to the angle <f) (= ~) 
with k — q (|k — q| ~ |k|). The incoming beam is polarized 
either perpendicular or parallel to scattering plane, and the 
scattered is a summation of both. The qy axis is taken to 
the same direction as the in-pane component of k; i.e., the 
illustrated qn is pointing the positive direction. 



3z 2 



component of both the minority and majority 



spin states at the energy ~ 0.3-E r below the Fermi level 
and a rich xy component of the majority spin just above 
the Fermi level illustrated in Ref. [7j. This coincides with 
the results here. Thus, we conclude that the 0.3E r mode 
in the AFM case simply arises from the interband tran- 
sition. 

In the spin-longitudinal mode, the flat excitation 
branch near 0.2E r shows zx — > xy character, and the dis- 
persive part of this branch near (§,0) possesses xy — > xy 
character as well. 

The excitation properties discussed above are expected 
to be observed by Fe L-edge RIXS experiments. In Fig. El 
we illustrate the geometry of our RIXS calculation. The 
momenta of incoming (k) and scattered (k') x-ray are 
perpendicular to each other (k _L k'), and these norms 
are almost the same k k' « 0.45267T in units of inverse 
lattice constant, where the lattice constant is assumed 
as 3.996 Ai^ The relation between the incident angle 9 
and the momentum transfer q = k — k' is determined 
as 9 = ^ + arcsin— with <f> (= ^), the angle between 
— k and k'. The incoming beam is either a polarized 
(perpendicular to the scattering plane) or ir polarized 
(perpendicular to the a polarization), and the scattered 
is a summation of both. 

The L-edge RIXS spectrum is calculated as a 
second-order perturbation involving the x-ray absorp- 
tion and emission processes. The transition operator 
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FIG. 4. (Color online) Fe Ls-edge RIXS spectra. Spectra 
for AFM case are plotted along (— f ,0)-(f ,0) in (a) and (e), 
(-f ,-§)-(! ,f) in (b) and (f), and (0 -f )-(0,f ) in (c) and 
(g); for the PM case, spectra are plotted along (— §,0)-(f-,0) 
in (d) and (h). Polarization is set as a for the upper panels 
and 7r for the lower. 



Dk for the electron system in the absorption process 
is given, within the dipole approximation, as ss 

E„,, M , CT ,k' c ^( £ ) 4' )(Ul aPk+k',J,i* + Hc > where P is the 

annihilation operator of the Fe 2p electron, and the 
coefficient is the dipole transition matrix cf^J (e) = 
(3d; n, a\e ■ r\2p;j,j z ), where j and j z represent the total 
angular momentum and its z component of the 2p elec- 
trons, respectively; e is the unit vector of the polarization 
of the incoming or outgoing x-ray. 

Assuming that the scattering occurs very fast, we em- 
ploy the fast-collision approximation; 20 i 22 which neglects 
the energy dispersion of the intermediate state. The spec- 
trum can be yielded from the calculated susceptibility 
and dipole transition matrix as 

Irixs(<1\\,u) oc -Im{ V VV 1 " 2 (q||,w) 

x [ E 4% ( £ i) c ii (^Ki (£o)] }, (3) 



where the subscript i (o) of e denotes incoming (outgoing) 
x ray; si and S2 are associated with (ax, a[) and (02, a' 2 ), 
respectively; and we here set j = 3/2, considering the Fe 
L3-edge absorption^ Due to the spin-orbit coupling of 
the Fe 2p electrons in the intermediate state, a spin flip 
can be achieved in the final state of the RIXS process 
— the dipole transition matrix elements between the op- 
posite spin states give finite contribution to the spectral 
intensity. This allows us to analyze the spin-transverse 
mode apart from the charge excitation. 

The calculated spectra for the AFM case are plotted in 
Figs.0Ja)-(c) and (e)-(g) along different directions of q|| 
for different polarizations. The fundamental structure of 
the spectra reflects the excitations discussed above, with 
the spectral intensity modified. What kind of excitation 
is enhanced depends on the polarization and q|| through 
the coupling of the spin and orbital degrees of freedom. 
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The magnon branch appears with a somewhat weak in- 
tensity, while the 0.3E r mode appears still clear. 

Note that the magnon observation in the present sys- 
tem is under the same condition as that in cuprates dis- 
cussed in Ref. In cuprates, the x 2 — y 2 spin can flip 
in the L-edge RIXS process because the local spin mo- 
ments lie in the Cu02 plane. In iron arsenides, the local 
magnetic moments also lie in the FeAs plane, where the 
dominant component, xy orbital, forms another linear 
combination of states with the z component of the or- 
bital angular momentum l z = ±2. Thus, the magnon in 
iron arsenides is also able to excite in the i-edgc RIXS 
process. 

We here describe the details of q|| dependence and po- 
larization dependence. The magnon branch is rather vis- 
ible along the diagonal path [Fig.^b) and (£)]: better in 
the +/— direction for the a/ir polarization. The 0.3E r 
mode tends to be more distinct for the negative . This 
mode along (— £,0)-(-|,0) shows a spin-transverse feature 
(spin-flip) for the a-polarization case, and shows spin- 
longitudinal and charge features (no spin-flip) for the tt- 
polarization case. 

In RIXS experiments using a twinned sample, the spec- 
tra will be obtained as a superposition of x and y direc- 
tions. In the spectra along (0,— ■|)-(0, : |), the basic struc- 
ture of the excitation spectra is similar to that along 
(— |,0)-(f ,0): magnon, 0.3E r mode, and 0.2E r mode. 
A difference appears near (— ^,0)/(0, — |), e.g., around 
0.2E r [see Figs. 2] (a) and (c)]: this anisotropy arises 
from the yz/zx orbital difference — these excitations near 
(— ^,0) show yz — > yz character. The spin character of 
0.3E r is also different. For example, the spin-longitudinal 
(charge) component contributing to the 0.3EV mode in 
the negative (positive) q|| makes the spin-transverse com- 
ponent no longer dominant in the spectrum for a po- 
larization. For 7r polarization, on the other hand, the 
spin-longitudinal component does not exist in the 0.3E r 
mode. To clarify the spin and orbital characteristics, 
RIXS experiments are required to be performed under 
the detwinned condition. 

The spectra for the PM case are plotted in Fig. EJd) 
and (h). There are no distinct excitation structures near 
0.3E r , in contrast to the AFM case. 3z 2 — r 2 — >• xy exci- 



tations lie near 0.2E r in the PM case; this orbital char- 
acter is consistent with the calculated density of states 
(not shown). The low-energy appearance is also differ- 
ent from the AFM case, spreading in the entire region. 
The q|| dependence of the spectral intensity appears in tt 
polarization. 

In summary, we have investigated the spin and or- 
bital characteristics of excitations in iron arsenides with 
a mean-field five-band Hamiltonian by RPA. In addition, 
we have calculated the L 3 -edge RIXS spectra with a 
fast-collision approximation. The spin character to be 
observed becomes selective by making use of the mo- 
mentum transfer dependence and the polarization depen- 
dence. We have discovered the spin and orbital charac- 
teristics expected to be observed for certain conditions 
as summarized below. 

We have reached a thorough understanding of the or- 
bital characteristics of the magnon excitation in the AFM 
state, revealed that the magnon excitation mainly in- 
volves the xy orbital, and concluded that the xy orbital 
plays a key role in the local magnetic moment in both 
the static and the dynamical spin features. We have 
confirmed that the magnon branch appears in L-edge 
RIXS spectra. Since it is weaker than orbital excita- 
tions in the RIXS spectra, it is useful to observe along 
(0,0)-(± : |,±f ), where +/— is for the <j/tt polarization. 

As for orbital excitations, a distinct branch has been 
found around 0.3E r in the AFM case, determined as a 
particle-hole excitation from 3z 2 — r 2 to xy. The differ- 
ence between the AFM and PM states should be observed 
as the absence of this mode. We propose an L-edge RIXS 
experiment under a detwinned condition to selectively 
observe the spin character of this mode. It will show a 
spin-flip feature for a polarization but no spin-flip for 7r 
polarization when observed along (— : |,0)-(0,0). 
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